Energy transmission in the forbidden bandgap of a nonlinear chain 
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A nonlinear chain driven by one end may propagate energy in the forbidden band gap by means 
of nonlinear modes. For harmonic driving at a given frequency, the process ocurs at a threshold 
amplitude by sudden large energy flow, that we call nonlinear supratransmission. The bifurcation 
of energy transmission is demonstrated numerically and experimentally on the chain of coupled 
pendula (sine-Gordon and nonlinear Klein-Gordon equations) and sustained by an extremely simple 
theory. 



Introduction. A nonlinear chain of oscillators has a 
few striking fundamental properties, the first of which, 
known as the Fermi-Pasta-Ulam recurence phenomenon 
[jl| is the spectral localization of energy: injected initially 
in one given eigenmode energy does not eventually dis- 
tributes amongst higher modes (as one would normally 
expect from a quasi linear approach). The FPU discov- 
ery has been at the origin of nonlinear studies all over the 
world and led to the birth of the soliton concept Q , soon 
followed by the creation of the inverse spectral tranform 
(1ST) and the concept of integrability |§ [Q §]. 

A nonlinear chain also shows up spatial localization of 
energy in the form of nonlinear coherent structures, the 
solitons [Q. This universal behavior is well understood 
in the concept of integrability: any localized bounded 
initial condition eventually evolves to a number of iso- 
lated solitons and a vanishing background of quasi-linear 
radiation. 

Nonlinear energy localization is now deeply studied in 
the context of discrete systems where intrinsic nonlinear 
modes (breathers for short) have been shown to be the 
fundamental basic objects [|, and have been experi- 
mentally observed in Josephson ladders |10| ], Moreover, 
as an effect of non-integrability, nonlinear modes do ex- 
change energy and the spatial localization acquires novel 



quite interesting features 11 



The behavior of a nonlinear medium submitted to 
boundary data, as opposed to initial conditions, is also 
of fundamental interest. One famous instance is the 
project of data transmission in optical fibers in nonlinear 
Kerr regime for which the basic model is the nonlinear 
Schrodinger equation [jT^. Another instance is the non- 
linear property of self-induced transparency (SIT) of a 
two-level system submitted to high-energy incident (res- 
onant) laser pulse JlJ] . In those two cases the concept 
of integrability has proved its efficiency as the physical 
boundary value problem maps to a well posed Cauchy 
(initial value) problem. 

There, the nonlinear coherent structures emerge from 
sufficiently energetic localized input pulse. Another fun- 
damental question is then the behavior of a medium sub- 
mitted to continuous wave radiation. Such problems have 
been studied so far mainly through peturbation of inte- 
grable models by external driving, see e.g. |l5| for a re- 
view and [ELm for interesting recent developpments. But 



the more basic problem of scattering of continuous wave 
onto nonlinear media has not been much studied. Pre- 
liminary numerical results have been obtained recently in 
p7| , where a linear monochromatic wave is scattered on 
a nonlinear medium. The nonlinearity allows then wave 
transmission under nonlinear modes generation, with a 
threshold which was not given a theoretical ground. 

Those problems share a common basic question: the 
response of a nonlinear medium to periodic boundary 
data. We demonstrate here the existence of a bifurca- 
tion of wave transmission within a forbidden band gap 
(FBG) in a nonlinear chain forced (periodically) at one 
end. The related brutal energy flow through the medium 
is illustrated on figure [l] where the energy E penetrating 
the medium is plotted as a function of the amplitude A 
of the boundary driving at a frequency in the FBG (see 
later for details). 




FIG. 1: 

for T = 



Plot of the energy E injected in the sine-Gordon chain 
140 as a function of the driving amplitude A. 



The mechanism of this bifurcation takes its origin in 
nonlinear modes generation by the periodic forcing and 
allows for energy injected at one end of a chain to pen- 
etrate a medium by an intrinsic nonlinear process. This 
is called nonlinear supratransmission and it is illustrated 
both by numerical simulations of sine-Gordon and of non- 
linear Klein-Gordon equations, and by experiments on a 
mechanical chain of pendula. 

We shall first prove that, for sine-Gordon, the mecha- 
nism for nonlinear modes generation follows a very simple 
rule providing an explicit formula for the bifurcation di- 
agram in the parameter space {A, tt} where A is the 
amplitude of the driving and fl its frequency. The rule 
states that energy penetrates the medium as soon as the 
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amplitude A of the harmonic driving at frequency fl ex- 
ceeds the maximum amplitude of the static breather of 
same frequency. Then the system adapts to the driving 
boundary by means e.g. of a 2-breather solution whose 
second component cannot be also static (a 2-breather so- 
lution of sine-Gordon is made of 2 breathers of different 
velocities) and hence propagates in the medium. The 
whole process eventually repeats itself. 

Nonlinear supratransmission holds also for the nonlin- 
ear Klein-Gordon chain obtained by Taylor expansion of 
sine-Gordon. Then this process does not relie on integra- 
bility and is expected to be generic as soon as the model 
possess a natural forbidden band gap. 

We discover then that harmonic phonon quenching en- 
hances nonlinear penetration : when the first significant 
harmonic generated by the nonlinearity (the third one in 
sine-Gordon) falls inside the FBG, the related phonons 
stick on the boundary and contribute to the driving al- 
lowing supratransmission at lower amplitude. This is an 
interesting property in view of experimental applications. 

Last the energy transmission by means of nonlinear 
modes generation is explored and the bifurcation is de- 
scribed by expressing the energy injected in the medium 
in terms of the driving amplitude, hence furnishing a 
striking view of the supratransmission process. 

Model. Consider the discrete sine-Gordon chain of 
coupled oscillators u n (f) (time is normalized to the eigen- 
frequency of the individual oscillator) 



c 2 (u n+1 



2u n + u n -i) +sinu n = 0, (1) 



on a semi-infinite line n > with given initial-boundary 
value problem, namely the data of the driving boundary 
uo(t), the initial positions u n (0), initial velocities w„(0) 
and boundary condition at the chain end. The linear 
dispersion relation u>(k) is given by 



1 + 2c 2 (l - cos k). 



(2) 



The chain will be submitted to external harmonic forc- 
ing uo(t) — Asm fit on a medium initially at rest. For 
a frequency £1 in the phonon band, quasi-linear waves 
are generated in the medium and, for large enough am- 
plitude, these waves will undergo Benjamin-Feir instabil- 
ity hence creating localized excitations. These nonlin- 
ear modes have a very important role in the large time 
asymptotic properties of a nonlinear system and are sus- 
pected to be responsible for turbulent-like behavior |L8j. 

We consider here a driving frequency in the FBG, 
namely Q < 1, for which the linear theory would lead 
to the evanescent wave Asm(£lt) exp[— An] with A given 
by 



A = arccosh ^1 



1 - fl 5 

2c 2 



(3) 



In the nonlinear case, in order to fit the boundary con- 
dition (ph, the medium adjusts a static breather 



Ub(n, t) = 4 arctan 



Ac sin(f2i) 



$1 cosh(A(n + no)) 



(4) 



which is an exact solution in the continuous limit only 
but works well for strong coupling (the fully discrete case 
c <C 1 will be considered separately) . 

Bifurcation process. Adjusting a static breather ac- 
tually means to adjust the value of the breather center 
—no such that the oscillation amplitude at the boundary 
n = matches the forcing amplitude. This works up to 
the maximum value A s of the breather amplitude real- 
ized for no = 0. Hence from ([|), the threshold A s reads 
as the following function of the frequency £1 



A, = 4 arctan 



(5) 



which has the accurate simplified continuous approxima- 
tion 4arctan[Vl — f2 2 /f2]. 

This bifurcation threshold is now checked on numeri- 
cal simulations of (|IJ) with the following initial-boundary 
conditions 



u (t) = Asm fit, u n (0) = 0, u n (0) = ASle~ 



(6) 



The initial velocities are those of an evanescent wave such 
as to partly avoid the shock wave generated by vanishing 
initial velocities (the same results, but time consuming, 
are obtained for vanishing initial velocities and a driving 
amplitude smoothly growing from the value to A). Fi- 
nally, an infinite medium is simulated by an absorbing 
boundary. The simulations are made with the dsolve 
routine in MAPLE with 10 5 maximum iterations. The 
absorbing end consists in adding a damping 'yii n in the 
model, with intensity 7 slowly varying from to 2 on the 
last 50 particules. 
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FIG. 2: Function U n (t) for n — 60 in the case fi = 0.90. 
The amplitude are A = 1.78 < A s for the first figure and 
A — 1.79 = A s for the second. 

Results. In order to generate a bifurcation diagram 
one has to compare between simulations where the non- 
linear supratransmission does or does not occur, as illus- 
trated on figures || where the motion of one particule of 
the chain is plotted for driving amplitudes just below and 
just above the threshold. The simulation is performed 
with 200 particles with a coupling c = 4. Each large os- 
cillation in the secong figure corresponds to a breather 
(constituted of a single hump oscillating and propagat- 
ing) passing by. Two of them are generated and cross the 
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site 60 at times 120 andl60. The small oscillation seen 
between the humps are the harmonic phonons, mainly of 
frequency 3fi. 

We may now proceed with a systematic exploration of 
the chain response. The result is presented on figure |3| 
obtained for 200 particles with a coupling c = 10 (some 
experiments have been actually made with smaller cou- 
pling and less number of points to shorten computation 
times) for a typical time of 200 (for frequencies close to 
the gap value 1, time had to be increased up to 500). The 
points on figure |^ are obtained with an absolute precision 
of 10~ 2 for the amplitude A. They are compared to the 
theoretical threshold expression (|j) (continuous curve). 




FIG. 3: Bifurcation diagramm in the (A, fi) plane. The curve 
is the graph of formula (p|). The points indicate the lowest 
amplitude for which nonlinear supratransmission starts. 

Harmonic phonon quenching. The figure |3| shows ex- 
cellent agreement to formula (||) at least in the region 
0.34 < O < 1. Discrepencies are seen to occur starting 
below 0.33 and 0.18. This results from the driving which, 
thanks to the nonlinearity, generates phonons at multiple 
frequencies (here third and fifth) . If these frequencies lie 
in the phonon band, the phonons move away from the 
boundary and have no effect on the forcing. If however 
they lie in the FBG, the related phonons do not propa- 
gate (which we call phonon quenching) and stick on the 
boundary where they add contribution to the driving. 

This effect should thus disappear if the phonons are 
eliminated by driving the boundary with the exact 
breather expression (Q) used to calculate u (t), u n (0) 
and u n (0). In that case we have checked that nonlinear 
supratransmission never occurs at an amplitude A < A s , 
while it occurs for very small deformation of the perfect 
breather. For instance, by using (1 + e)ttj,(0, t) as driving 
boundary, in the case fl = 0.30, supratransmission occurs 
for e = 6 10~ 4 , i.e. for a driving amplitude A — 5.0674 
instead of the threshold A s = 5.0644. 

Energy transmission. The nonlinear supratransmis- 
sion allows energy to flow through the medium and 
we compare here this energy flow below and above the 
threshold. The energy injected in the medium by the 
driving boundary is given at time T by 



In our case Uo(t) is the driving (0) and the chain is sup- 
posed infinite with u n (t) — > as n — > oo. Chosing for 
T an integer multiple of the period of excitation makes 
this energy to vanish identically in the linear case if the 
driving frequency falls in the FBG. 

In the nonlinear case, expression (Q) is computed nu- 
merically. For a driving frequency 0.9 and amplitudes 
running from 1.5 to 2.0, we obtain the figure |l| where the 
bifurcation is seen to occur for A = 1.80, the value pre- 
dicted by fromula ([s]). This simulation has been runned 
for frequencies in the range [0.2,0.9], with expected re- 
sults. 
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FIG. 4: Energy E injected in the Klein-Gordon nonlinear 
chain for T = 140 as a function of the driving amplitude A. 

Nonlinear Klein-Gordon. The approach stems from 
the existense of a breather solution of the model equation, 
allowing to determine the threshold amplitude. Then 
a fundamental question is the dependence of nonlinear 
supratransmission on some integrability property of the 
equation. To give an indication that this process is 
generic (with a stop gap), we have performed numeri- 
cal simulations of the following nonlinear Klein-Gordon 
chain: 



-c 2 (u„+i - 



- 2u n + u r . 



1 3 1 5 



0, (8) 



E = -C 



dt u (t)[ui(t) ~ Uo(t)} 



(7) 



the Taylor truncated expansion of sine-Gordon (the fifth 
order is kept to ensure a confining potential at large u n ). 

Then this system is solved with the boundary driving 
@ and the energy (0) is computed for the same param- 
eter values as for figure [|. The result is displayed on 
figure f|. Nonlinear supratransmission is seen to still oc- 
cur, though a nonlinear mode solution of the model does 
not exist. By scanning the frequency range in the gap, 
we have obtain that the process occurs down to f2 = 0.7 
and then disappears. 

Experiment. The phenomenon of nonlinear supra- 
transmission can be experimentaly realized on a mechan- 
ical pendula chain driven at one end by a periodic torque 
ifiTjf . The detail analysis of such experiments will be pub- 
lished later, but it is worth showing here a picture of 
a breather generated by the boundary driving at a fre- 
quency inside the FBG. The breather on picture || has 
been obtained with a chain of 30 pendula of angular 
eigenfrequency ojq — 15 Hz (upper value of the FBG) 
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by driving at (angular) frequency 12.7 Hz, which in the 
normalized units used here corresponds to SI = 0.85. The 
coupling constant has been measured to be c = 32. 




FIG. 5: Picture of a breather generated in a mechanical pen- 
dula chain driven at one end at a frequency in the forbidden 
band gap. 

Conclusion. A novel fundamental property of a non- 
linear medium has been unveiled, namely the capacity 
to transmit energy under irradiation in a forbidden band 
gap by means of nonlinear mode generation. Theoreti- 
cal construction of the bifurcation diagram is extremely 
simple when the one-breather solution is known and the 
numerical simulations fit strikingly well the theory. 



There remain of course many interesting open ques- 
tions and currently under study are the modulation of 
the driving signal, the effect of damping, viscosity, disor- 
der, external bias, discreteness, etc... Another essential 
point is the search of nonlinear supratransmission in oth- 
ers physical situations. For instance this result might 
provide understanding of the very mechanism of the gen- 
eration of gap solitons in photonic band gap materials 




Last but not least, as the sine-Gordon chain is a model 
for discrete Josephson transmission lines f22|| , we expect 
interesting applications in this field. There, the bound- 
ary driving is realized by microwave irradiation through 
a finline antenna [ p3| , and it actually corresponds to pre- 
scribing the derivative at the origin (Neuman condition). 
Preliminary numerical simulations has shown that non- 
linear supratransmission also works and that the bifurca- 
tion process obeys the same type of simple rule. Detailed 
results on this question wil be published later. 
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